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ABSTRACT 

We report the discovery of two additional planetary companions to WASP-41 and WASP-47. WASP-41 c is a planet of minimum 
mass 3.18 ± 0.20 Mj^p and eccentricity 0.29 ± 0.02, and it orbits in 421 ±2 days. WASP-47 c is a planet of minimum mass 1.24 ± 
0.22 Mjup and eccentricity 0.13 ± 0.10, and it orbits in 572 ± 7 days. Unlike most of the planetary systems that include a hot Jupiter, 
these two systems with a hot Jupiter have a long-period planet located at only ~1 au from their host star. WASP-41 is a rather young 
star known to be chromospherically active. To differentiate its magnetic cycle from the radial velocity effect induced by the second 
planet, we used the emission in the Htt line and find this indicator well suited to detecting the stellar activity pattern and the magnetic 
cycle. The analysis of the Rossiter-McLaughlin effect induced by WASP-41 b suggests that the planet could be misaligned, though 
an aligned orbit cannot be excluded. WASP-47 has recently been found to host two additional transiting super Earths. With such an 
unprecedented architecture, the WASP-47 system will be very important for understanding planetary migration. 

Key words, planetary systems - stars: individual: WASP-41 - stars: individual: WASP-47 - Techniques: photometric - Techniques: 
radial velocities - Techniques: spectroscopic 


1. Introduction 


In contrast to the large number of multiple planetary systems, 
stars with hot-Jupiter planets have long been thought to lack ad¬ 
ditional planets. Recent studies suggest that this statement may 
not be true (e.g. Knutson et al.|201^ . However the existence of 
those additional planets was revealed by nothing more than ra¬ 
dial velocity trends. Only seven outer planetary companions of 
close-in (a < 0.1 AU) giant planets have a full orbital period that 
has been observed. 

The multiple planetary system around v Andromedae ( |But-| 
iler et al.|1997| l was the first to be found with a hot Jupiter. This 
planet is surrounded by three more massive giant planets orbiting 
between 0.8 and 5 au (|Curiel et al.|201 1||. Ten years later, the dis¬ 
covery of the planetary system around HIP 14810 (Wright et al. 
2007 |2009|l revealed that a hot Jupiter can be found in a system 


with smaller planets on wider orbits. From 15 years of radial ve- 


* Using data collected at ESO’s La Silla Observatory, Chile: HARPS 
on the ESO 3.6m (Prog ID 087.C-0649 & 089.C-0151), the Swiss Euler 
Telescope, TRAPPIST, the L54-m Danish telescope (Prog CN2013A- 
159), and at the LCOGT’s Faulkes Telescope South. The data is publicly 
available at the CDS Strasbourg and on demand to the main author. 


locity surveys, [Feng et alT ( 2015| l recently constrained the orbits 
of the distant (P > 9 years) giant planets in orbit around the hot- 
Jupiter hosts HD 187123 and HD 217107. Among all systems 
with a close-in giant planet discovered by transit surveys, only 
HAT-P-13 dBakos et al.|2009ll, HA T-P-46 ( [Hartman et al.|20141 l, 
and Kepler-424 ( Endl et al.||2014[ ) have an additional planetary 
companion detected by Doppler surveys with a fully measured 
orbit. 


The Wide Angle Search for Planets (WASP, |Pollacco et al.| 
2006| l has discovered more than 100 hot Jupiters. While some 
of them clearly have a very distant companion (e.g. WASP-8, 
IQueloz et al.|2010l l, none of them were previously known to have 
a fully resolved orbit. 

The radial velocity follow-up observations of WASP initi¬ 
ated in 2008 with the fibre-fed echelle spectrograph CORALIE 
(mounted on the 1.2-m Euler Swiss telescope at La Silla) has 
identified new multiple systems including a hot Jupiter. We re¬ 
port here the detection of additional planets around WASP-41 
dMaxted etal.|2011| l and WASP-47 dHellier et al.|20T2| . These 
new planets are both located only ~lau of their parent star. 
When searching for long-period companions using the radial ve¬ 
locities technique, the effects of stellar activity should be care- 
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fully considered to avoid misdetections. A stellar magnetic cycle 
may mimic the signal of a planetary companion (e.g.|Lovis et al. 


|2011^ . WASP-41 is known to be active ( |Maxted et al.||2011| l, 
which makes it important to recognise the temporal structure of 
its magnetic cycle. Since the photospheric Ca (H & K) emission 
lines cannot be extracted from single CORALIE spectra (insuffi¬ 
cient signal-to-noise ratio), we alternatively use the emission in 
the Ha band to characterise activity. We find that this indicator 
is very well suited to WASP-41, though it may not be true in 
general ( |Cincunegui et al.|2007[ ). 

Section [^reports the observations, analysis, and results for 
WASP-47. Section similarly addresses the case of WASP-41 
with a description of the various additional observations, a study 
of the stellar activity, and an analysis that includes the Rossiter- 
McLaughlin effect followed by the results. In section after 
discussing the importance of considering magnetic cycles when 
searching for long-period planets, we compare our two multi¬ 
ple planetary systems, including a hot Jupiter, to the five previ¬ 
ously known ones. We discuss the biases induced by the observ¬ 
ing strategies of Doppler and transit surveys in finding multiple 
planetary systems including a hot Jupiter. Finally we check the 
transit probabilities of WASP-47c and WASP-41c. 


used in our analysis to characterise the transiting planet WASP- 
47b were 

- the transit depth defined as the planet/star area ratio (dF = 
(Rp/R*)^), where Rp and R* are the planetary and stellar ra¬ 
dius, respectively; 

- the transit impact parameter in the case of a circular orbit 
(b' = flb cos/p,b/R*), where Ub and /p_b are the semi-major 
axis and the inclination of the planetary orbit, respectively; 

- the transit duration T 14 , from first to last contact; 

- the orbital period Rb; 

- the time of inferior conjunction Ro.b; 

- the two parameters cos cob and yib sin Wb, where Cb is 
the eccentricity and cob the argument of periastron; 

- the parameter K 2 ^b = -^1 - where Kb is the radial 

velocity semi-amplitude. 

For the period Rb and time of inferior conjunction Ro,b, we im¬ 


During the process of revision of this paper, (Becker et al. 
2015| l reported the presence of two additional super-Farths tran¬ 
siting WASP-47 every ~0.8 and 9 days. This remarkable discov¬ 
ery brings strong constraints on the migration of this system and 
will be discussed in the last section. 


2. WASP-47 c 

2.1. Observations 

WASP-47 is a G9V star hosting a giant planet with a mass of 
1.14Mjup and a period of 4.16 days (|Hellier et al.||2012|l. The 
discovery paper includes 19 radial velocity measurements taken 
with CORALIE during two observing seasons. Ongoing moni¬ 
toring revealed a possible third body, and after 46 observations 
over 4.44 years, the orbit of a second planet was clear. After 
subtracting the inner planet (hot Jupiter) from the radial velocity 
data, the analysis of the residuals using a Lomb-Scargle peri- 
odogram shows a clear peak around 550 days (see Fig. [1]). The 
analysis of the residuals after subtracting the two-planet model 
shows no additional signal (see Fig. 0 . 

In November 2014, an upgrade of CORALIE was performed. 
The changes introduced an offset in the radial velocities. Since 
the orbit of the outer planet was not fully covered in phase, 
we collected six additional measurements from June to August 
2015. One measurement was done right after the upgrade just 
before the star became unobservable but it suffered from day¬ 
light contamination and had to be discarded. For the analysis, 
the 2015 data are considered as an independent radial velocity 
data set. 


posed Gaussian priors defined by the values derived in Hellier 
|et al.| ( |2012| l, since we did not include the WASP photometry 
in our analysis. Uniform priors were assumed for the probabil¬ 
ity distribution functions of all the other jump parameters. The 
choices of jump parameters were optimised to reduce the depen¬ 
dencies between the parameters and thus to minimise the corre¬ 
lations. We checked that using Ve cosw and yfe sinw as jump 
parameters translates into a uniform prior on e and does not in¬ 
fluence the results. 

We modelled the limb-darkening with a quadratic law, where 
instead of the coefficients u\ and M 2 , we used the combinations 
Cl = 2 xmi-i-M 2 andc 2 = mi- 2 xm 2 as jump parameters in order to 
minimise the correlation of the obtained uncertainties ( |Holman| 
|et al.||2006] l. We assumed normal priors on mi and M 2 with the 
values deduced from the tables by |Claret & Bloemen| ( |201 1} . As 
described in Gillon et al. ( 2012| l, we applied a correction factor 
CF — 1.32 to the error bars of the photometric data to account 
for the red noise. 

WASP-47 c does not have transit observations, the jump pa¬ 
rameters used in our analysis to define the second orbit were 

- the orbital period Rc; 

- the time of inferior conjunction Rq.c; 

- the two parameters y^coscUc and y^sinwc; 

- the parameter ^^ 2,0 = R^c “ 


1/3 


2.2. Analysis 

We combined the photometry data from EulerCAM (Gunn-r' fil¬ 
ter) published in Hellier et al. ( 2012| l with all the radial veloci¬ 
ties from CORALIE and used an adaptive Markov chain Monte- 
Carlo (MCMC) fitting scheme to derive the system parameters of 
the two-planet model. We did not include the WASP photometry 
in our analysis because it has poor precision compared to Eu¬ 
lerCAM and might suffer from spatial dilution from background 
stars. We used the most recent version of the code described in 
IGillon et al.|(|2012]) and references therein. The jump parameters 


We assumed uniform priors on all those jump parameters. 

The stellar density is derived at each step of the MCMC from 
Kepler’s third law and the values of the scale parameter Ob/R* 
and orbital period of WASP-47b ( Winn et al.||2010|. The stel- 
lar mass M* is obtained thanks to a calibration ( [Enoch et al.| 
|2010 Gillon et al.|201 l| l from well-constrained binary systems 
( [South worth |20 11 [ 1 . This empiric law is a function of Tee, P* 
and [Fe/H]. We used Gaussian prior distributions for Reft and 
[Fe/H], and the initial value for p*, based on the values derived 
by |Mortier et al. ( 2013| l. 


2.3. Results 

The results presented in Table have been obtained running 
five chains of 100000 accepted steps. The derived values are 
the median and 1 cr limits of the marginalised posterior distribu¬ 
tions, considering a 20% burn-in phase. The acceptance fraction 
is ~25% (~40% for the burn-in phase). We used the statistical 
test of Gelman & Rubin ( 1992| l to check the convergence of the 
chains and get a potential scale reduction factor (PSRF) smaller 
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Fig. 1. WASP-47 Lomb-Scargle periodogram of the residuals after sub¬ 
tracting Planet b. The dotted lines correspond to a false-alarm probabil¬ 
ity of 0.1%, 1%, and 10%. 


than 1.01 for all the jump parameters. The autocorrelations are 
shorter than 2 000 steps except for the semi-amplitude and the 
eccentricity of the outer planet. The lack of data around the ra¬ 
dial velocity minimum introduces a small degeneracy between 
these parameters. Their autocorrelations are ~5 000 steps long. 
As we did not include additional transit light curves, we did not 
improve the determination of the transit parameters compared to 
[Hellier et al. ( |2012| l. They are not given here for this reason. The 
best fit solution for the radial velocities is plotted in Fig. The 
values provided in Table [T are the median values of the poste¬ 
rior distributions and not the best fit values used in Fig. To 
verify that the median values represent a self-consistent set of 
parameters, we compare the BIC (Bayesian Information Crite¬ 
rion) obtained using the median values and the best fit values. 
We obtain compatible results. 

Including the data collected in 2015 allows to reduce the 
uncertainty on the eccentricity and semi-amplitude of the outer 
planet. All the other parameters remain unchanged compared to 
the results obtained with the first data set alone. The posterior 
distribution of the offset between the two data sets (ARV =18 
+ 11 m s“^) agrees with what has been observed on bright radial 
velocity standards. 

The semi-amplitude of WASP-47 c is barely twice as large 
as the radial velocity uncertainties, and the phase coverage is 
incomplete. To validate the detection, we compared the results 
obtained when running the same MCMC with a model including 
only one planet. With the one-planet model, the reduced on 
the initial set of radial velocities is 2.7 and the residual r.m.s. 
19.2 m s '. When fitting two planets, the corresponding reduced 
X^ is 1.2 and the residual r.m.s. 10.6 ms*. The two-planet model 
is clearly preferable. 


3. WASP-41 c 


3.1. Observations 

WASP-41 is a G8V star hosting a giant planet with a mass of 
0.9Mjup and a period of 3.05 days (|Maxted et al. 2011 1 . The 
discovery paper includes 22 radial velocity measurements taken 
within four months (excluding the first point). The short time 
span of the observations did not allow for detecting any addi¬ 
tional signal. We did a follow-up on the star and collected a 
total of 100 radial velocities with CORALIE over 4.46 years. 
The observations during the second season revealed an offset of 
~100ms * and triggered an intensive monitoring of the target. 
An additional planet on an eccentric orbit rapidly became appar¬ 
ent. Because of the significant eccentricity of the second planet 


Table 1. Median and l-tr limits of the posterior marginalised PDFs ob¬ 
tained for the WASP-47 system derived from our global MCMC analy¬ 
sis 


MCMC Jump parameters 

Star 


Peff [K] 

5576 ± 68 

[Fe/H] 

0.36 ± 0.05 

Cl 

1.17 ± 0.03 

^2 

-0.02± 0.02 

Planet b 


Pb[d] 

4.1591409 ± 0.0000072 

ro,b-2450000 [BJDtdb] 

5764.3463 ± 0.0002 

-^ib coswb 

0.04 ± 0.06 

sin Gjj, 

0.01 ± 0.09 

K 2 .b = K,.^l-elPl'^ [m s-' d'/^] 

225.7 ± 3.6 

Planet c 


Pc [d] 

572 ±7 

ro,c-2450000 [BJDtdb] 

5981 ± 16 

cos a>c 

-0.21 ± 0.24 

^e/ sinwc 

0.16 ±0.20 

P2.C = K, Vi^ [m s-‘ d‘/3] 

248 ± 43 

Derived parameters 

Star 


Wj 

0.462±0.014 

U2 

0.242±0.010 

Density p* [po] 

0.68 ± 0.06 

Surface gravity log g* [cgs] 

4.32 ± 0.05 

Mass M* [Mq] 

1.026 ±0.076 

Radius P* [Pq] 

1.15 ± 0.04 

Luminosity L* [Lq] 

1.14±0.11 

Planet b 


Mass Mp,b [Mjup] 

1.13 ±0.06 

RV amplitude Pb [m s“'] 

140 ±2 

Orbital eccentricity et, 

< 0.026 (at 2 a) 

Orbital semi-major axis [au] 

0.051 ± 0.001 

Density pp,b [pj] 

0.71 ± 0.08 

Surface gravity log gp,b [cgs] 

3.33 ± 0.03 

Equilibrium temperature [K] 

1275 ± 23 

Planet c 


Minimum mass Mp c sin ip ^ [Mj„p] 

1.24 ± 0.22 

RV amplitude [m s“*] 

30 ±6 

Orbital eccentricity 

0.13 ± 0.10 

Argument of periastron Wc [°] 

144 ± 53 

Orbital semi-major axis [au] 

1.36 ±0.04 

Equilibrium temperature [K] 

247 ± 5 

adR* 

255 ±7 


Notes. ^‘'Assuming zero albedo and efficient redistribution of energy. 


combined with a period close to one year, we had to wait until the 
fifth year of follow-up to plan observations at the periastron with 
a good sampling. After subtracting the inner planet (hot Jupiter) 
from the radial velocity data, the analysis of the residuals using a 
Lomb-Scargle periodogram shows a clear peak around 400 days 
(see Fig.[^. The additional peaks present at ~200 days and ~130 
days correspond to the two first harmonics P/2 and P/3. Such 
features are expected when applying frequency decomposition 
to eccentric signals, and WASP-4 Ic has an eccentricity e ~ 0.3. 
The analysis of the residuals, after subtracting the two planets 
from the radial velocities, shows no additional periodic signal 
(see Fig.|^. 

We used the HARPS spectrograph mounted on the 3.6-m 
telescope in La Silla to observe the Rossiter-McLaughlin effect 
induced by WASP-41 b during the transit on 3 Apr 2011. In addi¬ 
tion to the Rossiter-McLaughlin observations, 30 measurements 
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Fig. 2. WASP-47 CORALIE radial velocity data (blue/green dots: be¬ 
fore/after the upgrade) and best fit model (red line). Top: phase-folded 
on the period of the inner planet (outer planet subtracted). Bottom: 
phase-folded on the period of the outer planet (inner planet subtracted). 



Period (days] 


Fig. 3. WASP-41 Lomb-Scargle periodogram of the residuals after sub¬ 
tracting Planet b. The dotted lines correspond to false-alarm probabili¬ 
ties of 0.1%, 1%, and 10%. 


Table 2. Additional transit observations of WASP-4 lb and applied cor¬ 
rection factors. 


Instrument 

Filter 

date 

CF 

FTS 

PanSTARRS z 

12 Apr 2011 

1.8 



21 Mar 2011 

1.0 



02 Apr 2011 

2.1 

TRAPPIST 

I + z 

12 May 2011 

2.0 



09 Mar 2012 

1.5 



19 Apr 2013 

1.4 

Danish 

R 

19 Apr 2013 

2.2 



25 Apr 2013 W 

3.0 


Notes. Partial transit 


out of transit of WASP-41 were gathered with HARPS between 
Apr 2011 and Aug 2012 to confirm the second signal detected 
with CORALIE. 

More transit light curves of WASP-41 b have been obtained 
since the publication of the discovery paper: with the Faulkes 
Telescope South (FTS) at Siding Spring Observatory and with 
the TRAPPIST telescope ( |Jehin et al.|201 l| l and the Danish tele¬ 
scope both in La Silla. Details of these observations can be found 
in Table |2] 
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3.2. Study of stellar activity 


WASP-41 is known to be an active star ( |Maxted et al.|[20TT] l. 
In spite of the large amplitude of the second signal, we check 
that it is not due to a stellar magnetic cycle. For that purpose we 
searched for correlations between the radial velocities and spec¬ 
tral activity indicators. No variation is observed in the bisector 
span of the cross correlation function (Fig. [^bottom panel). The 
full width at half maximum (FWHM) shows some variability, 
which is suspiciously similar to the radial velocities’ variability 
(periodicity close to 400 days). The logR^jj. extracted from the 
HARPS spectra, based on the photospheric Ca II (H & K) lines, 
indicates a mean value of -4.483 + 0.036 characteristic of an ac¬ 
tive star. But no variation in the log is visible during the two 
years of HARPS data (Fig. [^. 

The emission in the Ha line is a rather underused activity in¬ 
dicator to trace photospheric stellar activity. In their study, |Cin-| 
cunegui et al. ( |2007| l concluded that the activity measured in 
the Ha was not equivalent to the indicator. They did not 
measure similar correlations between the two indicators in their 
whole sample. Some of their stars showed correlations, when 
others showed no or anti-correlations. Using solar observations, 
[Meunier & Delfosse] ( |2009[ ) demonstrated that the Ca II (H & 
K) and the Ha emissions have different behaviours. During an 
active period when the filaments reach saturation, they observe 
a positive correlation between the indicators. On the other hand, 
during moderate active periods, the filaments counteract the ef¬ 
fect of plages on the Ha measurements, leading to weak or neg¬ 
ative correlations. [Gomes da Silva et al.| ( |2014] l conclude that 
active stars with a mean logR^jj^ > -4.7 systematically show a 
positive correlation between the Ca II (H & K) and Ha flux. 

Since WASP-41 is a very active star with logR^jj. > -4.7, it 
is reasonable to use Ha as an activity indicator. The motivation 
to use this band is that it is located in the red part of the stellar 
spectrum where a high signal-to-noise ratio can be obtained from 
the CORALIE spectra. As suggested in jCincunegui et aL] ( |2007| l, 
we measured the Ha emission at 6562.808 A and the continuum 
centred at 6605 A and av eraged on a window 20 A wide. Unlike 


Cincunegui et al. (2007 1 we extracted the Ha from a 0.6 A wide 


band instead of 1.5 A. 

The Ha index reveals a time variation similar to the one ob¬ 
served on the FWHM (Fig. that can be modelled by a third- 
order trend in time. It is then reasonable to conclude that we are 
observing the magnetic cycle of WASP-41. After subtracting the 
third-order trend, we analysed the Lomb-Scargle periodogram 
of the residuals in the FWHM. A clear peak appears at exactly 
one year (see Fig. |^. The same one-year signal in the FWHM 
has been observed on the brighter stars also regularly observed 
with CORALIE, as well as with similar instruments (SOPHIE, 
HARPS). When we look at the Lomb-Scargle periodogram of 
the Ha index after subtracting the long-term trend, a peak ap¬ 
pears at ~18 days (see Fig.|^. Interestingly, this corresponds to 
the rotation period of the star derived from the WASP photome¬ 
try by jMaxted et al.| ( |201 l| l. 

Encouraged by the interesting result obtained using the 
Ha index on WASP-41, we looked at the similarly active star 


CoRoT-7 (with logR' = -4.62 > -4.7). Queloz et al. (2009 1 


note that CoRoT-7 shows strong variability in the FWHM, com¬ 
parable to what is observed on WASP-41. CoRoT-7 benefitted 
from simultaneous photometric and spectroscopic observations 
(see [Haywood et al.|2014 i. We extracted the Ha index from the 
HARPS spectra of CoRoT-7 in the same way as WASP-41. A 
very clear correlation is observed between the Ha index and the 
FWHM (see Fig.|7]l. After subtracting a third-order trend, very 
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Fig. 4. Magnetic cycle of WASP-41. Top: FWHM of the CCF from 
the CORALIE measurements. Bottom; Ho- index extracted from the 
CORALIE spectra. The red lines correspond to the fitted third-order 
trend. 
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Fig. 5. WASP-41 Lomb-Scargle periodogram of the FWHM after sub¬ 
tracting the third-order trend. A clear peak is present at 1 year. The 
dotted lines correspond to a false-alarm probability of 0.1%, 1%, and 
10%. 


0.30 
i 0.25 

a. 0.20 

T3 
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Fig. 6. WASP-41 Lomb-Scargle periodogram of the Ho- index after sub¬ 
tracting the third-order trend. A peak is present at ~ 18 days, correspond¬ 
ing to the photometric rotation period. The dotted lines correspond to a 
false-alarm probability of 0.1%, 1%, and 10%. 
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likely due to the magnetic cycle, a signal appears at ~23 days, 
corresponding to the photometric rotation period of the star. This 
similar result confirms that the Ha index is a good activity indi¬ 
cator for such active stars. 

The 400-day signal detected in the radial velocities is de¬ 
tected neither in the FWHM nor in the Ha. Its origin is most 
likely due to the presence of a second planetary companion. 


Fig. 7. Correlation between the Ha index and the FWHM from HARPS 
data of CoRoT-7. 


3.3. Analysis 

As we did with WASP-47, we simultaneously fitted the tran¬ 
sit light curves and the radial velocity measurements with an 
MCMC algorithm. To define the priors for the limb-darkening 
coefficients of the non-standard TRAPPIST I + z filter, we took 
the average of the values of the standard filters Ic and z', and 
the quadratic sums of their errors. We did not include the FTS 
photometry from the discovery paper in our analysis because it 
suffered from poor weather conditions. As suggested by Gillon| 
|et al.| ( |2012| l, we included a quadratic polynomial in time to 
model the transit light curves (except the fourth TRAPPIST and 
the first Danish transit observations). 

The observations with TRAPPIST on 21 Mar 2011 and 9 
Mar 2012 suffered from malfunctions of the autofocus. As a con¬ 
sequence, the FWHM of the images varied a lot. To correct for 
this effect, we added a third-order and a quadratic polynomial 
in FWHM, respectively, to model those two light curves. The 
TRAPPIST observation on 19 Apr 2013 experienced a meridian 
flip that we modelled including an offset for that light curve. The 
TRAPPIST light curve from the 2 Apr 2011 and the Danish one 
from the 25 Apr 2013 obviously exhibited stellar spot crossings. 
We visually identified the affected parts of the light curves and 
decreased their weights by increasing the corresponding errors. 
In general, we applied a correction factor to the error bars of the 
light curves to account for red noise (values listed in Table. |^. 

To improve the fit of the radial velocities and decrease the 
effects of stellar activity, we included a quadratic polynomial 
in the FWHM of the radial velocities from CORALIE. Then 
we quadratically added a "jitter" noise of 13.2ms“' to the er¬ 
ror bars to equal the mean error with the standard deviation of 
the best-fit model residuals. For the HARPS data (not includ¬ 
ing the RM sequence), we added a linear dependency in FWHM 
and a quadratic polynomial in bisector. We also added a "jitter" 
noise of 6.7 ms' to the radial velocity error bars. The baselines 
in FWHM and bisector were determined by comparing the BIC 
and the residual jitter obtained when running short chains with 
various models. 

We mode lled the Rossiter- McLaughlin effect following the 
equations of Gimenezj (2006 1 . We used Vv sin/* cosJ3, and 
Vv sin/* sinyS as jump parameters, where J3 is the projected an¬ 
gle between the stellar spin axis and the planet’s orbital axis, 
and V sin/* is the projected rotational velocity. The Rossiter- 
McLaughlin sequence is fitted as an independent dataset to mit¬ 
igate the effect of stellar activity. We did not need to add any 
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clT[h] 


Fig. 8. WASP-4 lb transit data (blue: FTS z'; green: TRAPPIST I + r, 
red: Danish R) and best fit model (black line). The TRAPPIST light 
curve was obtained by binning the data from five transit observations, 
and the Danish one from one full and one partial transit observations. 
Each point corresponds to a 7.2 min bin. 


"jitter" or dependency on other parameters to these data. Owing 
to the low transit impact parameter b^, there was a strong degen¬ 
eracy between p and v sin /* (see Triaud et al.||201 1 Anderson 


et al.|20lT]l. We imposed a normal prior on v sin /* to avoid un¬ 


physical solutions. Indeed from an MCMC with no prior on the 
V sin/* we obtainy 6 = -48°±29°, and v sin /* = 3.5^[ Q°kms“* 
with V sin /* getting higher than 100 km s ', clearly inconsistent 
with the spectral analysis. The value derived for p is strongly de¬ 
pendent on the value of the prior for v sin /* 

The priors used for Teff and [Fe/H] are defined with the val¬ 
ues from [Mortier et alT (2013|l. For the v sin/* we revised the 


value derived by Doyle et al. 
ters from [Mortier et al.| ( |2013 


( 2014| l, using the stellar parame- 
I to estimate the macroturbulence. 
Therefore we used v sin/* = 2.66 + 0.28kms * as a normal 
prior. Since we benefitted from transit observations that span 
more than three years, we did not use a prior on the period and 
time of inferior conjunction To b of the transiting planet. 


3.4. Results 

We present in Table the median and 1 cr limits of the 
marginalised posterior distributions after running five chains of 
100000 accepted steps (considering a 20% burn-in phase). The 
acceptance rate is ~25% (~40% for the burn-in phase). The con¬ 
vergence of the chains is validated by the test ofjGelman & Ru-j 
|bin| ( |1992'| ) (PSFR better than 1.01), and the autocoiTelations are 
shorter than 2 000 steps. Including the second planet in the model 
undoubtedly improves the fit. The r.m.s. of the CORALIE data 
after the two-planet fit is ten times smaller than after a one-planet 
fit. The best fit solution is plotted in Fig. for the transit light 
curves and in Fig.J^for the radial velocities. Again, the values 
provided in Table |^are not the ones used in Figs. and We 
checked that they represent a self-consistent set of parameters by 
comparing the BIC using both sets of values (best fit and median 
of the posterior distribution) and get compatible results. 

The best-fitting model suggests a misaligned planet with a 
projected orbital obliquity of p - -28° + 13°. However, the un¬ 
certainties are large (see Fig. [TOj i, and the data do not allow an 
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Table 3. Median and 1 cr limits of the posterior marginalised PDFs ob 
tained for the WASP-41 system derived from our global MCMC analy 
sis 


MCMC Jump parameters 

Star 

Teff [K] 

5545 ± 33 

[Fe/H] 

0.06 ± 0.02 

Cl.;' 

0.819 ± 0.009 

C2,;' 

-0.230± 0.006 

Cl./+; 

0.85 ± 0.04 

C2.;+; 

-0.19± 0.04 

Ci,R 

1.08 ±0.01 

C2,R 

Planet b 

-0.09± 0.01 

Planet/star area ratio (/Jp//?*)^ [%] 

1.87 ±0.01 

b' = flb cos ip.b//?* 

0.10 ±0.06 

ThM 

0.1099±0.0002 

Pb [d] 

3.0524040 ± 0.0000009 

ro,b-2450000 [BJDtdb] 

6014.9936 ± 0.0001 

yib coswb 

0.002± 0.048 

sintUb 

0.059 ± 0.074 

K 2 .b = Kbyll-elPl^^ [m s-‘ d'/^] 

199.7 ± 2.3 

Vv sin/* cosp 

1.42 ±0.10 

Vv sin/* sin/l 

-0.78!°i« 

Planet c 

Pc [d] 

421 ±2 

ro,e-2450000 [BJDtdb] 

6011 ±3 

COSWc 

0.53 ± 0.02 

sincuc 

-0.07± 0.06 

K2.C = K, [m s-' d‘/3] 

676 ± 21 

Derived parameters 

Star 


«1,;' 

0.282 ± 0.005 

U2,z' 

0.256 ± 0.002 

"l,/+z 

0.30 ± 0.02 

"2,/+; 

0.25 ± 0.02 


0.416 ± 0.006 

“2,R 

0.252 ± 0.004 

Density p* |po] 

1.41 ±0.05 

Surface gravity log g* [cgs] 

4.53 ± 0.02 

Mass M* [Mq] 

0.93 ± 0.07 

Radius R* [/Jq] 

0.87 ± 0.03 

Luminosity L* [Lq] 

0.64 ± 0.04 

Projected rotation velocity v sin /* 

Planet b 

[km s-‘] 2.64 ± 0.25 

Mass A/p,b [Mjcp] 

0.94±0.05 

Radius Rp,b [Rjup] 

1.18 ±0.03 

RV amplitude Rb [ms“'] 

138±2 

Orbital eccentricity Cb 

< 0.026 (at 2 cr) 

Orbital semi-major axis Ob [au] 

0.040 ± 0.001 

Orbital inclination ipb [°] 

89.4 ± 0.3 

Transit impact parameter />b 

0.10 ±0.06 

Projected orbital obliquity p [°] 

90+14 

^^-10 

Density Pp,b [pj] 

0.56 ± 0.03 

Surface gravity loggp.b [cgs] 

3.24 ± 0.01 

Equilibrium temperature [K] 

1244 ± 10 

Planet c 

Minimum mass Mp,c sin ip^c [Mj^p] 

3.18 ±0.20 

RV amplitude [m s“'] 

94 ±3 

Orbital eccentricity 

0.294 ± 0.024 

Argument of periastron a>c [°] 

353 ±6 

Orbital semi-major axis Oc [au] 

1.07 ±0.03 

Equilibrium temperature [K] 

241 ±2 

flc/R* 

265 ±3 


Notes. *** Assuming zero albedo and efficient redistribution of energy. 
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Fig. 9. WASP-41 radial velocity data (blue: CORALIE, green: HARPS; 
red: HARPS RM sequence) and best fit model (black line). Top: Phase- 
folded on the period of the inner planet (outer planet subtracted). Mid¬ 
dle: Zoom on the Rossiter-McLaughlin effect. Bottom: Phase-folded on 
the period of the outer planet (inner planet subtracted). 


3 . 7 r I-1-1-1-r 



1 9I- ■ ■_■ ■ ■ - 

-60 -40 -20 0 20 40 60 

/3[°] 

Fig. 10. Posterior distribution of fi and v sin/* for WASP-41 b on loga¬ 
rithmic scale. We impose a prior of 2.66 ± 0.28 km s“* on v sin/*. The 
contours correspond to the 1-, 2-, and 3-cr confidence limits for each 
parameter, and include 39.3%, 86.5%, and 98.9% of the MCMC steps, 
respectively. 


aligned orbit to be excluded (P(^| < 20°) ~ 24%). The rotation 
velocity expected from the photometric rotation period of the 
star (~ 2.44km s“') is compatible within the error bars with the 
spectroscopic v sin/* (2.64 + 0.25kms“*), suggesting that the 
star is equator-on. That reinforces the possibility of an aligned 
planet. 


4. Discussion 

We have shown that finding long-period planets with instruments 
such as CORALIE is still possible around active stars. Recog¬ 
nising the periodic structure of activity-related signals is essen¬ 
tial to avoiding false detections. While spots or plages induce 
signals at the rotation period of the star, the evolution of their 
coverage of the stellar surface is correlated with the magnetic 
cycle of the star. Such cycles have periods of several years and 
must be distinguished from long-period companions. |Lovis et al.| 
( |2011| l studied the effect of magnetic cycles on the radial veloc¬ 
ity measurements from the HARPS planet search sample. They 
derived empirical models of correlations between and ra¬ 
dial velocities, FWHM, contrast, and bisector spans, depending 
on the effective temperature Teif and metallicity [Fe/H] of the 
star. However, their sample only contains ’quiet’ stars because 
they are the preferred targets for planet searches. They consider 
only a handful of targets with a mean log/?[jj^ > -4.7, so their 
results cannot be applied to active stars. The CORALIE search 
for planets observes stars with a broader range of log/?J^j^ than 
the HARPS one. While the emission in the Ca II (H & K) lines 
is hard to extract from the CORALIE spectra, we have shown 
that the emission in the Ha, which is easy to measure, is a good 
indicator of activity for stars with log > -4.7. A systematic 
study of the correlation between the Ha emission and the radial 
velocities of active stars could help to estimate the influence of 
magnetic cycles on radial velocity measurements of active stars. 

The transit technique has led to the discovery of more than 
200 hot Jupiters, while only ~30 have been discovered by 
Doppler surveys. A total of nine multiple planetary systems that 
include a hot Jupiter (with a < 0.1 au) and a second planet with 
a fully probed orbit, are known (source: exoplanet.eu). Among 
them, four were discovered by radial velocity surveys and five 
by transit surveys. Those multiple systems represent ~13% of 
the known hot Jupiters found by radial velocities and only ~2% 
of those found by transit searches. The most obvious expla¬ 
nation for this difference is historical. Radial velocity surveys 
started discovering hot Jupiters well before transit surveys. The 
hot Jupiters discovered by radial velocities have been followed 
for more than ten years, allowing the detection of very wide 
companions (P > 9 years). The dedicated large transit surveys 
(mainly HATNet and SuperWASP) started to efficiently discover 
hot Jupiters about six years ago, allowing only the characteri¬ 
sation of relatively close companions. If we exclusively consider 
companions within 2 au, the rate of hot Jupiters with companions 
from the Doppler sample is divided by two. 

The difference in the rates of detected additional planets in 
systems with a hot Jupiter may also be explained by the dif¬ 
ferent observing strategies used in Doppler and transit surveys. 
The discovery of a hot Jupiter based on radial velocity data re¬ 
quires a large number of observations because the observer has 
no prior information on the planet. About typically ten radial ve¬ 
locity observations are carried out to confirm the existence of a 
hot Jupiter previously detected by transit. In that case the ob¬ 
server knows the period of the planet and can plan the radial ve¬ 
locity observations in order to measure the radial velocity orbit 
very efficiently. This means that, on average, hot Jupiters discov¬ 
ered by the Doppler technique benefit from a larger number of 
radial velocity measurements than those discovered by transits. 
Increasing the number of radial velocity measurements makes 
the detection of longer period planetary companions easier. 

In the beginning of the WASP follow-up with CORALIE, we 
started a systematic search for long-period companions around 
the WASP targets with a confirmed hot Jupiter. We have followed 
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Fig. 11. Known multiple planetary systems with fully probed orbits in¬ 
cluding a hot Jupiter (a < 0.1 au). Triangles: hot Jupiter found by radial 
velocity surveys. Circles: Hot Jupiters found by transit surveys. Dia¬ 
monds: systems presented in this paper. The black dashed lines corre¬ 
spond to a radial velocity semi-amplitude of 10 m s“', 30m s“*, and 
80m s“' assuming a circular orbit around a I-Mq star. Most of these 
planets were detected only by radial velocities, thus we are plotting the 
minimum mass M sin i instead of the true mass M. 
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more than 100 WASP host stars for durations of two to eight 
years, including about 90 targets followed for more than three 
years. The publication of this work is in preparation (Neveu- 
VanMalle et al. in prep). In this sample only two stars with a hot 
Jupiter have been found to host a second planet that has com¬ 
pleted at least one orbit. We still have a lack of multiple systems 
within 2 au compared to the hot Jupiters discovered by radial ve¬ 
locity surveys. 

Detectability thresholds can contribute to explain the ob¬ 
served discrepancy. Figure [TT| shows the multiple planetary sys¬ 
tems (with fully characterised orbits) that include a hot Jupiter, 
distributed in a plot of semi-major axis versus mass. The sym¬ 
bols located on the left-hand side of the plot (a < 0.1 au) are the 
hot Jupiters, and the symbols on the right-hand side (a > 0.1 AU) 
are the long-period planets. The dotted lines were drawn to vi¬ 
sualise which planets are part of the same system. We notice 
that except for WASP-47 c, only the long-period planets from 
the Doppler-survey systems have semi-amplitudes smaller than 
80m s'. Indeed, Doppler surveys target small samples of bright 
stars, while transit surveys target large samples of fainter stars. 
Besides being brighter, radial velocity targets are in general less 
active than transit targets. This means that a better precision can 
be obtained on the radial velocity measurements of stars from 
Doppler surveys. With a higher signal-to-noise ratio, detecting 
smaller signals is easier around brighter quiet stars. Smaller 
long-period planets might exist around stars with a known tran¬ 
siting hot Jupiter, but they are out of reach of current surveys. 

Knowing that WASP-47b and WASP-41b are transiting, it 
is natural to wonder if their outer planetary companions tran¬ 
sit as well. If we assume that the inclination of the outer planet 
is uniformly distributed, independently of the inclination of the 
hot Jupiter, the transit probability is ~0.4% for WASP-47c and 
~0.6% for WASP-41c (using [Kane & von Braun|2008| l. If we as¬ 
sume a coplanar system where the orbit of the outer planet can be 
uniformly distributed around 5° from the orbit of the hot Jupiter, 
we obtain a transit probability of ~6% in both cases. Detecting 
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the transit of the outer planet would bring strong evidence that 


the hot Jupiter had migrated inside the disk via gap opening (Lin 
|et al.|19^|Marzari & Nelson|2009l l. However, a non-detection 
would not be sufficient to infer a high mutual inclination between 
the orbits of the planets, which would be evidence of a migration 


Fabrycky & Tremaine 

2007 

[Nagasawa et al.|2008 Matsumuraj 

et al.|20101|Naoz et al. 

2011 

1. 


The next predicted inferior conjunctions for WASP-47c are 
2 Nov 2016 + 30 d and 24 May 2018 + 43 d, with an expected 
transit duration of ~18 h. A radial velocity follow-up is ongoing 
in order to improve the ephemeris. The next predicted inferior 
conjunction for WASP-41c (02 Nov 2016 + 7 d) cannot be ob¬ 
served by ground-based facilities because the star is close to the 
Sun. The two following ones are 12 Jan 2018 + 8d and 3 Mar 
2019 + 10 d with an expected transit duration of ~13 h. 

Recently, ([Becker et al. 2015|l have announced the detec¬ 


tion of two additional transiting planets around WASP-47. We 
checked if the nine-day planet could be detected in our set of 
radial velocities. Unfortunately, the expected amplitude of the 
signal is below the sensitivity of CORALIE for such faint stars. 
The geometry of this system favours disc migration theories be¬ 
cause scattering is unlikely to have kept the inner system copla¬ 
nar. Sanchis-Ojeda et al. ( 2015| l have recently announced a low 
stellar obliquity for WASP-47 that reinforces the disc-migration 
hypothesis. Recent results have shown that Kozai interactions do 
not affect multi-planetary systems. In this case the outer planet 
could have ‘protected’ the inner planets from secular interac¬ 
tions. This system brings fundamental constraints for planetary 
formation and migration theories. 
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Julian date - 2 400 000 


Fig. 1. CORALIE data of WASP-47 plotted in time. Blue/Green dots; before/after the upgrade. Top panel: Radial velocities superimposed with 
the best fit model. Second panel: Residuals of the radial velocities after subtracting the hot Jupiter superimposed with the model for the second 
planet. Third panel: Residuals of the radial velocities after subtracting the two planets. Bottom panel: Bisector spans. 
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Fig. 2. Radial velocity data of WASP-41 plotted in time. Blue: CORALIE; Green: HARPS; Cyan: RM sequence. Top panel: Radial velocities 
superimposed with the best fit model. Second panel: Residuals of the radial velocities after subtracting the hot Jupiter superimposed with the 
model for the second planet. Third panel: Residuals of the radial velocities after subtracting the two planets. Bottom panel: Bisector spans. 
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Table 1. Initial values and priors for WASP-47 


Parameter 

Value 

Prior 

M* [Mo] 

1.07 + 0.10 ' 

> 0 

Teff [K] 

5576 + 68 ' 

Gaussian 

[Fe/H] 

0.36 + 0.05 ' 

Gaussian 

Ui 

0.4584 + 0.0144 ^ 

Gaussian, mi + M 2 < 1 

U2 

0.2417 + 0.0086 ^ 

Gaussian, mi + M 2 < 1 

Pb[d] 

4.1591399 ±0.0000072 3 

Gaussian 

Tob-2450000 [BJDtdb] 

5764.34602 ± 0.00022 ^ 

Gaussian 

[%] 

1.051 +0.014 3 

uniform > 0 

ri4[d] 

0.14933 +0.00065 ^ 

uniform > 0 

K 

0.14 ±0.11 ^ 

uniform on b' - ab cos ip,b/J?*, 0 < Zttr ^ OblR* 

Kb [ms^'l 

136.0 +5 3 

uniform on K 2 ,b - Kb- e^ and > 0 

eb 

0 

uniform on cos Wb and sin Wb, if > 1 then eb = 0.999 

COb 

0 

uniform on yib cos c+b and sin a>b 

Pc [d] 

571 + 8 

uniform Od < P < lOOOyr 

7’o,c-2450000 [BJDtdb] 

5983 + 26 

uniform 

Kc [ms^^] 

34+ 10 

uniform on K 2 ,c = A’c -\/l ~ el Pq ^ and > 0 

ec 

0.2 ± 0.2 

uniform on coswc and sinwc, if > 1 then ec = 0.999 

Wc 

0 

uniform on cos tOc and sin tOc 


References. ' From 


Mortier et al. 


12013 I. ^ Interpolated from 


Claret & Bloemen 


(2011 1 . ^ From 


Hellier et al. 


| 2 ^. - 


See 


Kipping (2013 I. 


Table 2. Initial values and priors for WASP-41 


Parameter Value Prior 

M* [Mo] 0.90 + 0.07 1 >0~ 

T^s [K] 5546 + 33 ' 


[Fe/H] 

0.06 ± 0.02 ‘ 

Ml.z' 

0.2817 + 0.0044 2 

tt2,z' 

0.2559 + 0.0022 ^ 

tt\,I+z 

0.3090 ± 0.0340 ^ 

U2J+Z 

0.2510 + 0.0170 2 

U\,R 

0.4144 + 0.0065 2 

U2,R 

0.2520 ± 0.0037 2 

V sin 7* [km s“^] 

2.66 ± 0.28 2 

Pb[d] 

3.052401 ± 0.000004 

Fob-2450000 [BJDtdb] 

6014.991 ±0.001 4 

(R^/R.f [%] 

1.86 ±0.04 4 

Tidd] 

0.108 + 0.002 4 

btr 

0.40 + 0.154 

Kb [ms^*] 

135.0 + 84 

eb 

0 

a>b 

0 

jsn 

0 

Pc [d] 

418 ± 3 

ro,c-2450000 [BJDtdb] 

6006 ± 10 

Kc [ms^'] 

94 + 4 

ec 

0.3 ± 0.05 

a>c 

0 


Gaussian 

Gaussian 

Gaussian, + U 2 ^' < 1 ^ 

Gaussian, ui^' + U 2 ^’ < 1 ^ 

Gaussian, mij+j + M 2 ,/+z < 1 ^ 

Gaussian, mi_/+j + M 2 ,/+z < 1 ^ 

Gaussian, u\ r + M 2 ,« < 1 ^ 

Gaussian, u\ r + M 2 ,« < 1 ^ 

Gaussian 

uniform0d<P< 1000 yr 
uniform 
uniform > 0 
uniform > 0 

uniform on b' - a\, cos ip,b/f?*, 0 <bi^ < abIR* 
uniform on K 2 ,b = ^b Pb^ > 0 

uniform on cos tOb and -^/ib sin Wb, if > 1 then eb - 0.999 
uniform on yib cos a>b and yib sin Wb 
uniform on s/v sin/* cos (3 and Vv sin/* sinyS 
uniform0d<P< 1000 yr 
uniform 

uniform on \/l - e? P\^^ and > 0 

uniform on cos Wc and sin Wc, if > 1 then e^ - 0.999 
uniform on cos Wc and sin Wc 


Refere nces. 

(20111. 5 See 


From ] 

dortier et al. 

2013 

Kippinj 

;(2013|(. 


(2013 I. ^ Interpolated from 


Claret & Bloemen 


(2011 1 . ^ Updated from Doyle et al. 


(20141.''From 


Max ted et al. 
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Table 1. New radial velocities for WASP-47 


BJD-2450000 

(d) 

RV 

(kms“^) 

CfRV 

(kms-') 

FWHM 

(kms-i) 

Bis 

(kms-i) 

6205.690749 

-27.16568 

0.01258 

8.54665 

0.01824 

6213.632998 

-27.10977 

0.01866 

8.50854 

-0.12360 

6246.539542 

-27.03050 

0.01077 

8.52734 

-0.00957 

6269.531829 

-27.10498 

0.01370 

8.60976 

-0.06453 

6270.530518 

-26.95272 

0.01043 

8.56773 

-0.01388 

6437.858225 

-27.00067 

0.00972 

8.52917 

-0.02456 

6452.924978 

-27.00623 

0.00890 

8.56443 

-0.06375 

6459.924071 

-27.21081 

0.00784 

8.56435 

-0.03869 

6477.721354 

-27.05162 

0.01149 

8.59406 

-0.03366 

6486.844362 

-26.92702 

0.01094 

8.52159 

-0.02556 

6505.757749 

-27.20042 

0.00939 

8.55957 

-0.06514 

6512.745209 

-26.99712 

0.01142 

8.49187 

-0.06260 

6514.752564 

-27.12983 

0.01053 

8.50852 

-0.04074 

6515.749341 

-26.95270 

0.01151 

8.53025 

-0.02725 

6530.629824 

-27.21972 

0.01850 

8.54674 

-0.02802 

6534.789870 

-27.22405 

0.01955 

8.57303 

-0.04356 

6559.713758 

-27.24359 

0.01907 

8.51742 

-0.05402 

6590.650676 

-26.97269 

0.01437 

8.51798 

-0.06681 

6616.574051 

-27.05365 

0.01745 

8.54661 

-0.04370 

6782.899610 

-26.99561 

0.01233 

8.53542 

0.00179 

6808.876282 

-27.19380 

0.01140 

8.54211 

-0.04406 

6831.876641 

-26.96141 

0.00732 

8.57188 

-0.06365 

6859.804354 

-27.16125 

0.00982 

8.48966 

-0.02312 

6885.631966 

-26.97427 

0.01222 

8.57285 

-0.00047 

6915.753320 

-26.94052 

0.01116 

8.51941 

-0.00544 

6950.553823 

-27.19106 

0.00775 

8.48781 

-0.04470 

CORALIE upgrade 




7175.831586 

-27.17962 

0.01397 

8.47192 

-0.05504 

7185.803516 

-26.94267 

0.01691 

8.53709 

-0.07024 

7200.842220 

-27.19041 

0.01242 

8.48536 

-0.04705 

7229.724940 

-27.21735 

0.02664 

8.49685 

-0.06773 

7253.631858 

-27.15060 

0.01096 

8.50363 

-0.06696 

7258.670969 

-27.24158 

0.01503 

8.50061 

-0.06009 



A&A-W41c_W47c_corr, 



5650 5700 5750 5800 5850 5900 5950 6000 6050 6100 6150 
JD - 2450000.0 [days] 


Fig. 3. WASP-41 extracted from the HARPS spectra. 
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Table 2. New radial velocities for WASP-41 


BJD-2450000 

(d) 

RV 

(kms-') 

O^RV 

(kms-1) 

FWHM 

(kms-') 

Bis log/^HK criog«;,K 

(kms-') 

CORALIE 





5577.833651 

3.53952 

0.01033 

8.81258 

-0.04113 

5591.760168 

3.30884 

0.00989 

8.80484 

-0.01806 

5595.802071 

3.50882 

0.01022 

8.81885 

0.00150 

5597.776083 

3.22930 

0.00773 

8.80704 

-0.02692 

5607.721786 

3.41489 

0.00873 

8.75296 

-0.01255 

5621.850786 

3.20236 

0.00929 

8.87451 

-0.00206 

5623.709905 

3.45711 

0.01054 

8.81959 

0.00608 

5624.829879 

3.20966 

0.00976 

8.80052 

-0.00711 

5625.826430 

3.34101 

0.00878 

8.78968 

-0.03361 

5629.892800 

3.45361 

0.00973 

8.86781 

-0.01814 

5635.815355 

3.50088 

0.01104 

8.82368 

-0.07613 

5646.754274 

3.22713 

0.00777 

8.79920 

-0.01367 

5647.630652 

3.46375 

0.00793 

8.83207 

-0.00998 

5649.657082 

3.19961 

0.00800 

8.82716 

-0.01678 

5675.784327 

3.38673 

0.01022 

8.78597 

0.01963 

5676.605991 

3.17071 

0.00771 

8.72380 

-0.00364 

5677.773655 

3.29207 

0.00895 

8.72341 

-0.04367 

5680.715317 

3.30325 

0.01022 

8.73045 

0.00763 

5691.683825 

3.25986 

0.01085 

8.78840 

-0.03964 

5705.592134 

3.43592 

0.01549 

8.76269 

-0.00381 

5716.577101 

3.14687 

0.00798 

8.68982 

-0.00490 

5764.502263 

3.32023 

0.01004 

8.60551 

-0.00142 

5769.525384 

3.33383 

0.00953 

8.73791 

0.03592 

5952.825523 

3.54587 

0.00779 

8.74638 

0.00973 

5953.844642 

3.46188 

0.01269 

8.69680 

0.00014 

5954.786428 

3.34226 

0.00838 

8.68532 

-0.00038 

5955.794145 

3.57514 

0.01921 

8.71114 

-0.01929 

5958.776062 

3.54863 

0.00813 

8.79174 

-0.00790 

5961.787850 

3.50127 

0.00799 

8.77488 

-0.00942 

5962.724308 

3.56518 

0.00845 

8.74688 

-0.02158 

5963.826866 

3.33818 

0.00795 

8.77194 

-0.02822 

5974.729812 

3.56877 

0.00869 

8.66496 

-0.03499 

5976.866742 

3.45179 

0.00929 

8.76283 

-0.01263 

5990.888836 

3.30386 

0.00775 

8.72585 

0.00909 

5991.827941 

3.31816 

0.00824 

8.71079 

-0.03771 

5992.810751 

3.56087 

0.00741 

8.74076 

0.00721 

6002.863880 

3.35757 

0.01043 

8.75194 

-0.00926 

6003.827288 

3.26345 

0.01026 

8.74074 

0.00006 

6004.617666 

3.47529 

0.01352 

8.74425 

-0.00198 

6031.771885 

3.34243 

0.01269 

8.65108 

0.00612 

6032.514592 

3.49173 

0.01228 

8.67336 

-0.01604 

6033.583037 

3.28047 

0.00955 

8.66582 

-0.01872 

6067.660548 

3.16248 

0.00874 

8.65557 

0.01180 

6068.614165 

3.36070 

0.01015 

8.65640 

-0.01275 

6113.472992 

3.14763 

0.01296 

8.66391 

-0.04187 

6115.524431 

3.34283 

0.00784 

8.62192 

-0.02416 

6116.552476 

3.16081 

0.00982 

8.61649 

0.01291 

6117.563578 

3.37099 

0.00935 

8.67662 

-0.02545 

6140.504929 

3.17964 

0.01699 

8.66987 

-0.09129 

6150.483360 

3.19956 

0.01465 

8.59768 

-0.02625 

6151.480248 

3.42741 

0.01413 

8.62581 

-0.00604 

6272.851056 

3.30816 

0.01512 

8.68548 

-0.01227 

6285.856902 

3.46249 

0.01426 

8.67411 

0.00109 

6311.817319 

3.26497 

0.00849 

8.70165 

-0.05757 

6314.770425 

3.26354 

0.00906 

8.69916 

-0.00225 

6335.760287 

3.38387 

0.01025 

8.69464 

-0.03972 

6344.893416 

3.37439 

0.01489 

8.67605 

-0.01501 

6362.834685 

3.50514 

0.01487 

8.69525 

-0.06191 

6364.806955 

3.54831 

0.01140 

8.71556 

-0.00035 
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Table 2. Continued. 


BJD-2450000 

(d) 

RV 

(kms-') 

O^RV 

(kms-1) 

FWHM 

(kms^*) 

Bis 

(kms^') 

log/^HK 

CTlogRJjK 

6365.700329 

3.53718 

0.01081 

8.74924 

-0.02150 



6372.694569 

3.30787 

0.01324 

8.71766 

-0.00158 



6401.508276 

3.52006 

0.00950 

8.71947 

-0.02537 



6419.727309 

3.46173 

0.00981 

8.62232 

0.01914 



6427.640513 

3.25900 

0.01451 

8.64472 

-0.02081 



6441.643745 

3.48042 

0.01110 

8.64098 

-0.01737 



6450.560548 

3.46452 

0.00870 

8.65585 

-0.02942 



6484.570287 

3.42376 

0.01376 

8.57077 

-0.01220 



6496.467235 

3.44219 

0.00958 

8.62008 

-0.01283 



6506.482624 

3.27505 

0.01018 

8.59371 

-0.00833 



6514.470015 

3.38330 

0.00963 

8.64153 

-0.02157 



6646.830447 

3.29643 

0.00838 

8.70805 

-0.02489 



6685.792726 

3.46022 

0.00978 

8.73929 

-0.03507 



6699.747165 

3.22170 

0.00926 

8.76697 

-0.03725 



6713.750758 

3.42864 

0.00862 

8.70401 

-0.03290 



6723.823895 

3.25117 

0.00909 

8.65259 

-0.01733 



6738.590042 

3.36799 

0.00905 

8.64089 

-0.04485 



6748.823614 

3.36232 

0.00913 

8.71895 

-0.04756 



6830.613737 

3.26410 

0.01257 

8.63168 

-0.05510 



HARPS 







5655.551700 

3.22759 

0.00363 

7.59140 

-0.00330 

-4.4430 

0.0172 

5655.853857 

3.23851 

0.00633 

7.58522 

0.00779 

-4.4359 

0.0332 

5656.536367 

3.39812 

0.00321 

7.61857 

0.02087 

-4.4596 

0.0154 

5656.697723 

3.43224 

0.00265 

7.57389 

0.00204 

-4.4619 

0.0116 

5656.862517 

3.45792 

0.00341 

7.59920 

0.02669 

-4.4865 

0.0194 

5681.502346 

3.47325 

0.00470 

7.54826 

-0.01900 

-4.4496 

0.0228 

5681.747792 

3.46028 

0.00331 

7.52105 

-0.02056 

-4.5356 

0.0190 

5682.508177 

3.27667 

0.00278 

7.54712 

-0.01070 

-4.4821 

0.0126 

5682.746435 

3.22961 

0.00357 

7.55328 

-0.00814 

-4.4868 

0.0178 

5683.607603 

3.29092 

0.00271 

7.53416 

-0.00548 

-4.4662 

0.0116 

5684.631995 

3.48142 

0.00302 

7.56635 

-0.01153 

-4.4767 

0.0136 

5685.587048 

3.27317 

0.00249 

7.57172 

0.00787 

-4.4737 

0.0104 

5711.539738 

3.39743 

0.00428 

7.55737 

0.00110 

-4.4846 

0.0208 

5712.569788 

3.38172 

0.00379 

7.55858 

-0.00532 

-4.5048 

0.0191 

5716.557913 

3.16014 

0.00336 

7.47622 

0.01511 

-4.6466 

0.0252 

5728.497874 

3.20106 

0.00352 

7.48383 

-0.00680 

-4.5084 

0.0169 

5729.498999 

3.27894 

0.00418 

7.50356 

-0.00211 

-4.4838 

0.0202 

5749.508251 

3.28400 

0.00444 

7.47420 

-0.01793 

-4.5368 

0.0250 

5753.559500 

3.19748 

0.01026 

7.50218 

-0.05764 

-4.4422 

0.0560 

5754.523986 

3.43646 

0.00761 

7.50424 

-0.02168 

-4.4780 

0.0417 

5776.464239 

3.41784 

0.00841 

7.49327 

0.00614 

-4.4815 

0.0489 

6020.679485 

3.49349 

0.00572 

7.45429 

-0.01218 

-4.5778 

0.0433 

6021.672381 

3.25952 

0.01203 

7.46629 

-0.05230 

-4.4977 

0.0799 

6029.682164 

3.47547 

0.00966 

7.45214 

-0.03355 

-4.5089 

0.0648 

6058.586834 

3.17781 

0.01144 

7.38666 

-0.03210 

-4.4681 

0.0662 

6059.587695 

3.43300 

0.01062 

7.46955 

0.01536 

-4.5694 

0.1073 

6092.564596 

3.24052 

0.00991 

7.47168 

-0.02810 

-5.0597 

0.2935 

6093.534543 

3.45736 

0.00745 

7.45766 

-0.00441 

-4.5485 

0.0567 

6094.550435 

3.26252 

0.01267 

7.42991 

-0.03793 

-4.5425 

0.1083 

6151.471858 

3.43600 

0.01277 

7.41152 

0.00263 

-4.5729 

0.1182 

HARPS RM sequence 






5654.541742 

3.44100 

0.00408 

7.58857 

-0.01057 

-4.4695 

0.0211 

5654.627495 

3.42446 

0.00309 

7.56772 

-0.00340 

-4.4528 

0.0134 

5654.694950 

3.40718 

0.00521 

7.56499 

-0.03342 

-4.4779 

0.0331 

5654.701224 

3.40040 

0.00512 

7.55081 

0.00192 

-4.5322 

0.0367 

5654.707323 

3.39846 

0.00446 

7.57235 

-0.02608 

-4.4858 

0.0272 

5654.713365 

3.39405 

0.00544 

7.54723 

0.00574 

-4.5201 

0.0379 

5654.719580 

3.39300 

0.00525 

7.58325 

-0.02580 

-4.5016 

0.0384 

5654.725854 

3.39530 

0.00487 

7.57183 

-0.00835 

-4.4374 

0.0310 

5654.731722 

3.39696 

0.00480 

7.57224 

0.01508 

-4.4085 

0.0283 
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Table 2. Continued. 


BJD-2450000 

(d) 

RV 

(kms-') 

O^RV 

(kms-1) 

FWHM 

(kms^*) 

Bis 

(kms^') 

log/^HK 

CTlogRJjK 

5654.738169 

3.38889 

0.00486 

7.58768 

-0.00323 

-4.5086 

0.0356 

5654.744152 

3.38091 

0.00449 

7.57718 

-0.03317 

-4.5097 

0.0331 

5654.750368 

3.38292 

0.00460 

7.62232 

0.01020 

-4.5330 

0.0359 

5654.756525 

3.38866 

0.00483 

7.58194 

-0.01348 

-4.4997 

0.0358 

5654.762741 

3.39363 

0.00506 

7.57987 

-0.01364 

-4.4967 

0.0373 

5654.768667 

3.40373 

0.00531 

7.55134 

0.00864 

-4.4188 

0.0320 

5654.774824 

3.41298 

0.00579 

7.56770 

-0.01484 

-4.4786 

0.0416 

5654.781317 

3.40804 

0.00595 

7.55392 

-0.01020 

-4.4952 

0.0432 

5654.787243 

3.38972 

0.00548 

7.58143 

-0.01048 

-4.5308 

0.0446 

5654.793586 

3.38582 

0.00517 

7.57592 

-0.03576 

-4.4906 

0.0375 

5654.799339 

3.38918 

0.00622 

7.58407 

-0.01323 

-4.4909 

0.0470 

5654.805785 

3.36643 

0.00549 

7.60594 

-0.00493 

-4.5604 

0.0485 

5654.811885 

3.36490 

0.00604 

7.61466 

0.00560 

-4.4852 

0.0462 

5654.818043 

3.35075 

0.00696 

7.61693 

-0.02194 

-4.4819 

0.0523 

5654.824385 

3.34658 

0.00659 

7.62860 

-0.00536 

-4.5539 

0.0568 

5654.835242 

3.33034 

0.00638 

7.60329 

-0.00396 

-4.5984 

0.0655 

5654.841399 

3.32018 

0.00596 

7.56792 

0.01640 

-4.6577 

0.0683 

5654.847557 

3.31924 

0.00596 

7.60258 

-0.03724 

-4.5498 

0.0518 

5654.853541 

3.33239 

0.00702 

7.58480 

0.02337 

-4.6217 

0.0726 

5654.859930 

3.34743 

0.00778 

7.56688 

-0.02332 

-4.4810 

0.0624 

5654.866029 

3.34820 

0.00757 

7.61977 

0.01059 

-4.6210 

0.0795 

5654.872071 

3.34226 

0.00772 

7.60082 

-0.03158 

-4.5566 

0.0707 

5654.878402 

3.33363 

0.00903 

7.59843 

0.00049 

-4.7133 

0.1221 

5654.884502 

3.33904 

0.00873 

7.57028 

0.00809 

-4.5821 

0.0908 

5654.890659 

3.34222 

0.00860 

7.59899 

-0.03549 

-4.3253 

0.0494 

5654.896829 

3.34804 

0.00859 

7.57287 

0.01621 

-4.5036 

0.0757 



